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ABSTRACT 

In this work we extend the study on the mass distribution of the spiral galaxy NGC 5278, 
performing ID and 2D (GALFIT) bulge-disk decomposition to determine which components 
constitute the baryonic mass in this galaxy. Our analysis does not detect any bulge, instead we 
find a bright source, probably related with the central AGN, and an exponential disk. We fix the 
stellar disk contribution to the rotation curve (RC) with broad band photometric observations and 
population synthesis models, to obtain the 2D mass distribution of the stellar disk. In particular, 
for NGC 5278, we find that the typical assumption of considering the mass-to-luminosity ratio 
(MjV) of the disk as constant along the galactocentric radius is not valid. We also extract a 
baryonic RC from the mass profile, to determine the inability of this baryonic RC (also taking into 
account ± 30% errors in the disk mass), to fit the entire RC. We perform the RC decomposition 
of NGC 5278 considering the determined baryonic RC and four types of dark matter (DM) halo: 
Hernquist, Burkcrt, Navarro, Frenk, & White and Einasto. Our results determine that Hernquist 
DM halo better models our observed RC in the case of disk mass Md = 5.6 x 10 10 M and also 
with less 30% disk mass. In the case of more 30% disk mass the cored Einasto (n < 4) DM halo 
is the best fitting model. 

Subject headings: Galaxies: interactions, (Cosmology:) dark matter, Galaxies: kinematics and dynamics, 
Methods: data analysis, Techniques: interferometric, Techniques: imaging spectroscopy 



1. Introduction. 

NGC 5278 is a spiral galaxy participating in 
an interaction process with the companion NGC 
5279. Both galaxies constitute the isolated galaxy 
pair KPG 390 that belongs to the Karachent- 
sev catalog of isolated pair s of galaxies in the 
northern hemisphere (KPG) ()Karachentsev 19721 ) 
and it is classified as a M 51 galaxy type by 
Klimanov fc Reshetnikovl (|200ll ). The pairs of 
galaxies in the KPG catalog are in the first stage 
of interactions and in particular KPG 390 being 
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considered M 51 galaxy type is not in a merger 
status. NG C 5278 and NGC 5279 are Markarian 
galaxies and iKeel et al.l ( 19851 ) have classified the 
two galaxies as intermediate type between Seyferts 
and liners, this means that the center of NGC 5278 
is the host of an Active Galactic Nucleus (AGN). 
Another important information we obtain from 
the literature is that the morphological classifica- 
tion of NGC 5278 is not unique: Sb (HyperledcQ 
(|Paturel et alJl2003l) 1. SO D (Simbad) and SA(s)b? 
pec (NED). 

We have already performed a detailed study 
of the interac t ion p rocess in a previous work 
(|Repetto et al.l l201oh . In our earlier study we 



derive velocity fields and RCs of both members of 
KPG 390. The hypothesis to build the RC from 



1 http: / /leda. univ-lyonl.fr 
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our Ha velocity fields is that in the relation of the 
line of sight velocity the terms of radial velocity 
and vertical velocity are small enough with respect 
to the terms of rotational velocity. In the case of 
NGC 5278 we have already encountered that con- 
dition inside an angular sector of 20° around the 
kinematic major axis of this galaxy (the residuals 
are less than 60 km s _1 within that angular sec- 
tor of 20°) and we derive the RC of NGC 5278 
within this angular sector to remove the effect 
of non-circular motions. Inside that sector, the 
RC of NGC 5278 really traces the gravitational 
potential required to decompose its RC and the 
corresponding parameters (V S y S , P. A. and inclina- 
tion) are free from any uncertainty derived from 
possible regions of non-circular motions in NGC 
5278. 

In general the RCs of spiral galaxies are important 
tools used to infer the contribution to the gravi- 
tational potential of the different components of 
a galaxy: stel lar disk, bulge and bar, dust, gas, 
and D M halo (|van der Kruitlll989t ISofue fc Rubin 
200lh . As a matter of fact, the existence of very 
massive DM halos was inferred from the flat- 
ness of the RCs at large galactocentric distances 



( Rubin et al. 1980). The kinematic data used to 



derive the RCs could come from absorption lines 
of the stars in the disk of the galaxy, from emis- 
sion lines of the ionized gas in HII regions of the 
disk, or from neutral gas in the disk. 
The method to study the DM halo distribution 
takes the RC, constructed from the kinematic ob- 
servations of a spiral galaxy, and expresses it as 
V 2 (R), where V 2 (R) is the square of the rota- 
tional velocity derived from the kinematics and 
R is the galactocentric distance. Then it fits the 
contribution to V 2 (R) of the assumed different 
components such as: a stellar (exponential) thin 
disk, a stellar bulge, the HI gas, and a DM halo 
with different possible radial density profiles. A 
non linear least squares minimization is used to 
select which disk and halo density distributions 
fit better the observed kinematics. Several impor- 
tant assumptions are made in using this method. 
Let us discuss the ones concerning the stellar disk 
contribution to the kinematics. One of the ma- 
jor problems arises because from the photometry 
we obtain only the surface brightness distribution, 
not the surface mass density distribution. This im- 
plies that we have to assume a certain M/L ratio 



of the stellar disk to convert surface brightness to 
mass density. In our earlier work we accomplished 
the RC decomposition of NGC 5278, considering 
a disk component without performing bulge-disk 
decomposition, and setting free the M/L ratio of 
the disk (within acceptable values), because we 
do not have any constraint of that value. In that 
former study we performed the RC decomposi- 
tion considering th ree DM halos: Hernquist halo 
(IHernquistl Il990al). Pseudo-I sothermal Halo (P- 
ISO) (|van Albada et al.lll985l) and Navarro, Frenk 



White Halo (NFW) (jNavarro et al.l I1996D . In 
iRepetto" et all (l2010l ) we were unable to discrimi- 
nate between the three DM halos employed to de- 
compose the RC of NGC 5278. The three halos re- 
produce the RC of NGC 5278 well enough with dif- 
ferent disk masses and scale lengths and the values 
obtained for the halo masses and radii are in ac- 
cordance with the typical values re ported by cos- 
mological numerical simulations ([Navarro et al 
19961) . 



In all the previous studies on galaxy mass determi- 
nation through RC decomposition, the Mj L ratio 
of the disk was considered constant along the ra- 
dius and taken as a free paramet er. One of the first 
appro aches was proposed by Ivan Albada et al" 



(|l985l ) who varied the disk M/L from its mini- 
mum value (zero, i.e., no stars in the disk) to its 
maximum value (the one that fitted the inner RC). 
This approach is k nown as "the maximum-disk 
fit". For instance, iBlais-Ouellette et al.l (jl999l ) 
study the Sd galaxy NGC 5585, considering the 
M/L ratio of the stellar disk as a free parame- 
ter, in the process of mass decomposition of the 
RC of this galaxy. The authors obtained that 
the P-ISO DM halo succeeded in reproducing the 
observed RC of NGC 5585, nonetheless the re- 
sulting M/L ratio of the stellar disk that better 
fits the different components had unrealistic low 
values (M/L ~ 0.3 — 0.8). Since then, much 



work (Ide Blok et all [20011: ISwaters et al I 12003 ; 



Spano et al. l2008t IChemin et al.l 2011 ) has been 



done on the mass distribution of galaxies through 
RC decomposition, always considering the M/L 
ratio as a free parameter, leaving the principal 
recipe (M/L ^const) unchanged; nevertheless, 
using this simple re cipe, the M/L ratio values are 



often unphysical. iFuentes-Carrera et al. ( 2007 ) 



have tried to fix the M/L of the disk with more 
reasonable values obtained from the disk color in- 
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dexes, comparing them with those of stars and 
adopting the M/L of those stars as stellar disk 
M/L. The authors obtained as best fit DM halo 
the P-ISO DM profile, however they consider that 
the disk is composed of stars of only one spectral 
type; instead a detailed stellar population synthe- 
sis study is necessary, to better constrain the disk 
M/L. 

The motivation of the present study is to ex- 
tend the RC decomposition of NGC 5278, this 
time with new constraints on the baryonic mass 
of this galaxy, derived from additional observa- 
tions concerning surface photometry in several 
band s and using stellar population synthesis stud- 
ies ( Zibetti et al.l l2009') to determine the baryonic 
M/L of NGC 5278. From the broad band sur- 
face photometry analysis and the stellar popu- 
lation synthesis we obtain the disk mass profile 
of NGC 5278. From that mass profile we derive 
the baryonic RC of NGC 5278 by means of the 
cumulative integration of the surface brightness 
along the line of sight, to establish the inadequacy 
of this baryonic RC and also of the baryonic RC 
more or less 30% disk mass to account for the en- 
tire mass in NGC 5278. In the mass modeling of 
NGC 5278 we also explore a disk mass variation 
of more or less 30% conside ring four types of DM 
halos: Hernq uist, Burkert (lBurkertl ll995h. NFW 
and Einasto (|Einastolll965l ). In this analysis we 
find that the RC of NGC 5278 can be reproduced 
well enough by Hernquist DM halo with a disk 
mass Md — 5.6 x 10 10 M Q and also in the case of 
less 30% disk mass. The cored Einasto (n < 4) 
DM halo produced the best fit to the RC of NGC 
5278 with more 30% disk mass. 
The article is divided as follows: In the introduc- 
tion we present the problem of fitting several com- 
ponents to galaxy RCs derived from kinematical 
observations. In section [2] we accomplish bulge- 
disk decomposition to establish the baryonic mass 
component of NGC 5278, showing that there is 
no bulge contribution. Section [3] addresses the 
method to obtain disk mass distributions from 
certain color indexes and population synthesis 
models. This section also gives the application 
of this method to obtain the stellar disk mass dis- 
tribution of NGC 5278. In Section 2] we describe 
the procedure to derive the baryonic RC of NGC 
5278 from the mass profile obtained in Section [3] 
In Section [5j once we have determined the true 



disk mass distribution, we use the baryonic disk 
of NGC 5278 to fit only the DM halo distribution 
in the observed RC. In Section |6] we consider the 
maximum disk solution for NGC 5278, where we 
utilize the baryonic disk determined in section 3] to 
fit the majority of internal points of the observed 
RC. In Section [7] a discussion of our results and 
their implications is presented and the conclusions 
are outlined in section |8] 

2. ID and 2D bulge-disk decomposition of 
NGC 5278. 

As we already mentioned in the Introduction 
the morphological classification of NGC 5278 is 
uncertain and in the Hyperleda database it is clas- 
sified as Sb, so we want to explore the possibil- 
ity of the presence of a bulge component in NGC 
5278. In the Introduction we have also reminded 
that NGC 5278 shows an AGN in its center, fact 
that could create more confusion in the morpho- 
logical classification. For these reasons we per- 
form ID bulge-disk decomposition of NGC 5278 
considering several photometric bands in three 
ground-based telescopes (DSS (B), SDSS (u, g, i), 
2MASS (J, K)) and in the HST (F300W, F814W). 
For all telescopes and photometric bands consid- 
ered we obtain the same result, in this section we 
only show the surface brightness profile fit to the 
2MASS K band and Hubble F814W. We made use 
of IRAlU task ellipse to obtain the surface bright- 
ness profiles along the major and minor axes of 
NGC 5278, a nd of IRAF task n fit Id to fit a de- 
Vaucouleurs (|de Va ucouleurs 1948) bulge and an 
exponential disk (lFreemanlll970l) accor ding to the 
procedure devised bv iKormend vl (|l977h . We make 
a PyRAF-based script to run both IRAF tasks for 
all the photometric bands considered. In this sec- 
tion we only show the surface brightness profile 
fit along the major axis because the profiles along 
major and minor axes are identical but the profile 
along major axis is more extended. The results are 
displayed in the Figure [T] and it is clear from the 
two panels that the de Vaucouleurs profile do not 
fit the inner part of the surface brightness profile 
i.e., the plausible bulge component of this galaxy. 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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Fig. 1. — Surface brightness profiles fit along major axis for the Hubble F814W and 2MASS K band. 
The "cdisk" nomenclature in the figures is referred to the cutoff exponential disk expressed by the relation 
J(r) = la cxp [— (ctr + /3 3 r~ 3 )], where /3 has length dimension s and I(r) reduces to an exponential disk for 
r ^ (3. This cutoff disk profile, according to iKormendvl (|1977| ). is necessary to better fit the underlying disk 
component. The zoom in the upper part of the figures shows how the deVaucoleurs fit to the expected bulge 
of NGC 5278 does not succeed. 
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Kormendv fc Kennicuttl (|2004l ) demonstrated 
that Sb and earlier-type galaxies contain prin- 
cipally classical bulges (deVaucouleurs), whereas 
Sbc galaxies have more often non-classical bulges 
(Sersic index between 1 and 2) and Sc and later 
type galaxies seem never to have classical bulges. 
These facts are very important and given the in- 
congruity in the morphological classification of 
NGC 5278 we are obliged to use a Sersic profile 
(|Sersid Il968l ) to fit the expected bulge compo- 
nent in NGC 5278. For this purpose we per- 
form 2D bulge-disk d e comp osition using GAL- 
FIlll (|Peng et al.ll2002L l2010h . The principal rea- 
son to utilize GALFIT is because that software 
makes a 2D decomposition of a galaxy image al- 
lowing to employ several components such as ex- 
ponential disk, Sersic, deVaucouleurs, Gaussian, 
King profile, Moffat, PSF and other components 
for bars, spiral arms and Fourier modes. The 
presence of a central AGN, i.e., a point source, 
and the fact that GALFIT cannot resolve sources 
with radius less than 0.5 pixels (corresponding to 
0.025 kpc for NGC 5278) and axis ratio less than 
0.1, forced us to perform the fit with the F814W 
Hubble image of NGC 5278 due to the high res- 
olution (0.0996" - 0.005 kpc for NGC 5278) of 
Hubble images and the fact that GALFIT needs 
high resolution images to work appropriately. The 
motivation to choose the F814W filter, which is 
equivalent to the Johnson I band, is substantially 
due to the lack of NICMOS infrared images for 
NGC 5278 that could better demonstrate the ex- 
istence of a classical bulge in NGC 5278. In fact, in 
the case of this galaxy there are only images from 
the WFPC2 in the three filters F814W, F255W 
and F300W. NGC 5278 is practically invisible 
in the filter F255W that is in the mid-UV, in 
the filter F300W there is considerable emission of 
this galaxy but F300W is in the near-UV, so the 
more suitable choice for our purpose (considering 
a classical bulge made up of old stars) is the filter 
F814W. We begin to fit a PSF function (fwhm 5.5 
pixels « 0.27 kpc) to eliminate the strong central 
compact peak and we succeeded to remove most 
of the central source (fwhm 7.4 pixels 0.37 kpc) 
proceeding then to fit an exponential disk compo- 
nent to the residuals of the first fit (85 counts). 
After removing other relics of the AGN we man- 



age to get residuals of almost 23 counts, similar to 
the counts of the spiral arms. In fact, the brighter 
spot of this last residual map is represented by an 
HII region in the periphery of the spiral arms (23 
counts). In these last residuals, the zone around 
the extracted central point source has 20 counts. 
Thus, the bulge-disk decomposition analysis in 
NGC 5278 do not detect any bulge component, 
instead we find a bright, almost punctual, central 
spot in the center, probably associated to the cen- 
tral AGN. The results of our GALFIT analysis are 
shown in Figure O 



3 http: / / users . obs . carnegiescience.edu / peng/work / galfit / galfit .html 
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Fig. 2.— GALFIT 2D decomposition of the NGC 5278 F814W Hubble image. In the residual image the 
spiral arms mean level is about 18 counts with a maximum of 23 counts, in this figure there is clearly visible 
the dust lane, more evident that in the original data where the disk emission of NGC 5278 partially hides 
this feature. In the same image the counts level around the central PSF is about 20. 
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3. Determining the disk M/L ratio of 
NGC 5278 from stellar population syn- 
thesis models. 

In order to determine the mass profile and to- 
tal mass estimate for NGC 5278 in an appropriate 
way, resolved stellar mass maps are needed. Struc- 
tural parameters may change between luminosity 
and stellar mass studies because the mass to lu- 
minosity ratio, M/L = (M* / M )/(£*/ Lq ) , is 
no t a constant throug hout the disk. The method 
of Zibetti et al.l (2009) is capable of reconstruct- 
ing resolved stellar mass maps of galaxies from 
multi-band optical/NIR imaging. The technique 
is based on a Monte Carlo library of 50000 stel- 
lar population synthesi s (SPS) models supported 
in the 2007 version of iBruzual fe Charlotl (|2003h 
code. These mode l s inclu de a new prescription by 
Marigo fe Girardil (|2007l ) for the thermally pul- 
sating asymptotic giant branch (TP-AGB) evo- 
lution of intermediate and low mass stars. The 
library covers a wide range of star formation his- 
tories, metallicities, and dust contents. The pre- 
scription to treat dust attenuati on is based on 
the tw o-component dust model of ICharlot fe Fall 
(2000). " n this prescription, the emission from 
stars younger than 10 7 yr, which is the time 
needed for the dissipation of dense molecular 
clouds, is more attenuated than that from older 
stars. 

We take advant age of the novel technique of 
Zibetti et al. ( 20091 ). and construct a spatially re- 
solved map of stellar mass surface density for 
NGC 5278. For this purpose we use the q-ban d 
and i-band data from SDSS DR8 (lAiharall201ll) 
and t he i^ s -band data from 2MASS (jSkrutskie 
19971 120061) . The procedure is as follows. The 
Monte Carlo spectral library was obtained from 
da Cunha etafl (|2008h . The MAGPHYS (Multi- 
wavelength Analysis of Galaxy P hysical Proper- 
ties) cod^ ( da Cunha et al.ll2008l) was used to ob- 
tain lookup tables of (i — K s ) color, (g — i) color, 
and M/L ratio at a redshift of z = 0.025 (see 
figure El). Models were binned in 0.05 magnitude 
bins, and the median value for each bin was com- 
puted. We compare the (i — K s ) and (g — i) colors 
obtained from the photometry with the values ob- 
tained from the binned models on a pixel-by-pixel 
basis. The resolved mass map was obtained by 



multiplying the K s luminosity image by the M/L 
ratio map. 

For all calculations we adopt the cosmolog- 
ical parameters Hq = 70.0 km sec -1 Mpc" 1 , 



a, 



0.3, and f2 



vacuum 



or Hi 



0.7. The lu- 



minosity distance adopted from NED was (m — 
M) = 35.24 mag, or 112 Mpc. The absolute 
mag nitude of the Sun in the K s -hand was taken 
from lBlanton fe Roweisl (|2007j) . No corrections for 
Galactic ex tinction were made since E(B — V) = 



1998). 



0.009 mag (jSchlegel et al. 

The resolved stellar mass map was depro- 
.jected adopting th e projection parameters of 
Repetto et al.l (|2010l) . The resulting mass to light 
(in the i^-band) ratio and the mass (azimuthally 
averaged) profiles are shown in figure 21 

4. The necessity of DM to decompose the 
RC of NGC 5278. 

From the mass profile of figure [5] (left panel) we 
derive the baryonic part of the disk of NGC 5278, 
considering the relation M(R) = 2ir f Q r E(R)RdR 
to compute the surface mass density inside the ra- 
dial extensions of the mass profile of figure [5] (33 
arcsecs w 17.3 kpc), then we calculate the circu- 
lar velocity with the relation V 2 (R) = GM(R)/R. 
We also analyze a variation of the baryonic part 
of the disk of NGC 5278 of more or less 30% disk 
mass, due to a 30% intrinsic uncertainties of the 
Zibetti et al. I (<2009l) method. The result is shown 
in figure [51 From this figure we observe that if 
we only consider the baryonic disk of NGC 5278 
there is needed at least 70% more mass to recover 
the RC of this galaxy. The fact that the whole 
RC of NGC 5278 requires the DM component to 
explain the mass distribution of this galaxy is not 
obvious from the fit we performed in our previous 
work, because in that earlier study we randomly 
varied the M/L ratio of the baryonic disk of NGC 
5278. In that former study we fitted jointly the 
free disk M/L component with the DM compo- 
nent and we obtained different values for the disk 
mass and scale radius for the different DM ha- 
los (P-ISO, Hernquist halo ( Hernquistl 1990a ) and 
NFW) we have analyzed (jRepetto et al. I I2010I) . 
That procedure determines the true parameters 



*http:/ /www. iap.fr/magphys/magphys/MAGPHYS. html 



5 The biggest uncertainty comes from propagating the error 
in the assumed distance to the object. 
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Fig. 3. — Effective M/L in i^-band as a function of (i — K s ) and (g — i) colors at a redshift of z = 0.025. 
Left panel: the median logarithmic M/L for models binned 0.05 x 0.05 mag 2 in the color-color space. Right 
panel: the rms of log M/L. 




Fig. 4 . — Left. Azimuthally averaged radial profile of M/L in JC s -band for NGC 5278, adopting I Zibetti et al 



(2009) method based on g, i, and K s photometry. Right: Radial mass profiles of NGC 5278 plotted in 
logarithmic scale, with the KPG 390 image superimposed to denote the regions where the profiles are taken 
in the galaxy disk of NGC 5278. The radial profiles were shifted vertically by the logarithmic amount of 10 9 . 
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Fig. 5. — Left: Mean logarithmic mass profile extracte d using the method outlined in section |3l Right: 
Ha observed rotation curve of NGC 5278, 
exponential function. 



obtained in iRepetto et al. (2010), preliminary fitted with an 




Fig. 6. — Baryonic disk of NGC 5278 more and less 30% disk mass. In this panel is also shown the total RC 
(disk+DM) of NGC 5278 for comparison. 
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of the baryonic component of NGC 5278, never- 
theless, the actual baryonic mass of NGC 5278 
could have been overestimated or underestimated, 
as well as the DM component of that galaxy. In 
the present analysis we avoid this problem con- 
straining the disk parameters by means of broad 
band surface photometry additional observations 
and stellar population synthesis studies so that we 
only set free the DM parameters. 

5. RC decomposition of NGC 5278 fitting 
only the DM component after subtrac- 
tion of the determined baryonic disk. 

The RC decomposition was performed using the 
module mpfil[^| (Levemberg-Marquardt method). 
We extract the DM RC of NGC 5278 subtract- 
ing the baryonic RC of that galaxy displayed in 
figure U from the observed RC of NGC 5278, also 
taking into account a disk mass variation of more 
and less 30% to consider a 30% uncertainty in 
the determination of the disk mass. The corre- 
sponding total baryonic mass estimates derived 
from the baryonic RC of NGC 5278 are respec- 
tively 5.6xl0 10 M© (disk mass), 7.3 xlO 10 M© 
(disk mass more 30%) and 3.9xl0 10 M© (disk 
mass less 30%). 

We proceed to explore the DM halo distribution 
by fitting the DM RC of NGC 5278. We consider 
this time four hal o models: Her nquist, Burkert, 
NFW and Einasto (jEinastolll965h ). In the present 
work we add two new halo profiles, the Burkert 
DM halo (|Burkertlll995h and the Einasto DM halo 
with three free parameters (|Einastol[l965l) . The 
reason to utilize Burkert DM halo instead of P- 
ISO is dictated by the fact that Burkert DM halo 
is more general than P-ISO and behaves like P- 
ISO (i.e. it presents a core) at smaller radii while 
converges to a NFW at larger radii. The expres- 
sion of the sp atial density o f the Burkert halo is 
the following (|Burkertlll995l ): 



Pb{t) 



(r + r )(r 2 + r§) : 



(1) 



where po and ro are free parameters that rep- 
resent the central DM density and the scale ra- 



6 We use the python Kmpfit vers ion from the Kapteyn pack- 
age l|Terlouw fe Vogclaar 2 012Tl t hat uses the C implemen- 
tation of mpfit Markwardt (2009). The corresponding web- 
site is http://purl.com/net/mpfit 



dius. The choice to include the Einasto DM halo 
is motivated by relatively recent A-CDM dissipa- 
tionless simulations of gravit ational clustering and 



tg ana 
_ l2004 

20101 : iTissera et al1l2010l) . This DM halo den- 



isolated spherical collaps e ([Navarro et al 



sity profile has already been used by other au- 
thors to fit the DM distribution when baryon s 
are included (lOnorbe et al.ll2007t[Gao et alJl2008h . 



Navarro et al.1 (|2004L l201fih proposed a model with 
functional form p (r) oc exp(— Ar a ), where a is a 
parameter describing the degre e of curvature of 
the profile. iMerritt et al. (2005) pointed out that 
this is the same relation previously proposed by 
Einasto to describe the sp atial density distribution 
of galaxies (|Einastolll965l) . It is important to note 
that this DM density profile introduces more free 
parameters so that, in principle, it could allow a 



better fit to the RCs. For example. iChemin et al 



(|2011l ) fit 20 RCs of the THINGS survey of isolated 
galaxies, with the Einasto, NFW and P-ISO pro- 
file and the result of their analysis is that Einasto 
index of less than 4 (cored Einasto) produces a 
better fit with respect to P-ISO in about 60% of 
the RCs considered, while Einasto index greater 
than 4 (cuspy Einasto) produces a better fit with 
respect to NFW in about 80% of the RCs. For 
all these reasons we use in this work the Einasto 
DM halo as well. IMerritt et al] (|2006l) provide the 
following formulation for the spatial density of the 
Einasto DM halo: 



PE (r) = p e exp {~d n [(r/r e )^ n - 1]}, (2) 

where p e is the density at the radius r e that de- 
fines a volume containing half of the total mass, 
the term d n is a function of n defined according to 
the approximation d n = 3n — 1/3 + 0.0079/n for 
n > 0.5. 

During the fitting procedure we leave free the DM 
halo radius and the DM halo mass (and also the 
Einasto index). In order to accomplish the fit- 
ting process, we required to start with initial DM 
halo guesses Mh and an, where Mh is the DM 
halo mass and ah is the DM halo radius. We 
vary these halos parameters within the intervals 
M h G[10 9 -10 14 ] M , a h £[0.1-40.0] kpc. For the 
Einasto DM density profile the range of variation 
of the halo shape parameter is n s[0. 6-10.0]. The 
choice of these intervals is twofold: the necessity 
to explore DM halo parameters domain as wide as 
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possible and the obligation to avoid possible re- 
gions where the parameters values could lose their 
physical meaning. We use these ranges for the 
DM halo masses and the DM halo radii also in the 
cases of more or less 30% disk mass. From the 
fitting process we obtain that the best fit to the 
DM RC is the Hernquist halo profile with a disk 
mass M d = 5.6 x 10 10 M© and less 30% of this 
disk mass value. For more 30% disk mass the best 
fitting model to the DM RC of NGC 5278 is the 
cored (n < 4) Einasto DM halo , nevertheless the 
value of the Einasto index, of our analysis, is out- 
side of the interval (5 < n S i m < 7) predic t ed b 



the numerical simulations of Navarro et al. (200 



2010). The fits to the RC for all disk masses and 



the four types of DM halos are shown in figures [3 
EH El H2 The values of the halo masses 
and halo radii are displayed in tables [TJ [3l [2] We 
can see from the masses values in these tables that 
the different mass models, in our study, have sim- 
ilar halo masses within a factor 4, and the con- 
centrations for the NFW halo and Einasto halo 
are in accordance with the last simula tion results 
(|Klvpin et al.ll201lL IPrada et al.ll2012h . From the 
three figures it is clear that the Burkert DM halo 
is not able to reproduce our observed RC, for the 
disk masses considered, and the final halo radius of 
3 kpc indicates that it is equivalent to P-ISO DM 
halo, because as mentioned above the Burkert DM 
halo has a flat core at smaller radii. 
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Fig. 7. — Rotation curve decomposition of NGC 5278 with disk mass Md = 5.6 x 10 




Fig. 8. — Rotation curve decomposition of NGC 5278 with disk mass Md = 5.6 x 10 
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Fig. 10. — Rotation curve decomposition of NGC 5278 with more 30% disk mass Md = 7.3 x 10 10 . 
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Table 1: Rotation curve decomposition with 
M d =5.6xl0 10 M Q . 



Principal Results 


Hernquist 


Burkcrt 


NFW 


Einasto b 


Halo Mass (M©) 


8.7X10 11 


5.5X10 11 


9.9X10 11 


1.6xl0 12 


Halo Radius (kpc) 


17.9 


3.0 


9.0 


18.6 


2a 
Ar 


0.99 


6.71 


1.28 


1.13 


DoP 


10 


10 


10 


9 


Additional Results 


Initial R vir (kpc) 


Final R vir (kpc) 


Initial c V i r 


1 mal c v i r 


NFW Halo 


150.0 


167.8 


10.2 


9.9 


Einasto Halo 


133.8 


198.1 


10.5 


9.6 


Reduced \ 2 an <3 degrees of freedom. 








'Final Einasto index n=3.6. 








Table 2: Rotation 


curve decomposition with 






M d =7.3xl0 10 M© (more 30% disk mass) 








Principal Results 


Hernquist 


Burkcrt 


NFW 


Einasto b 


Halo Mass (M©) 


1.3xl0 12 


3.3X10 11 


5.4xlO n 


l.lxlO 12 


Halo Radius (kpc) 


30.1 


3.0 


o o 
9.U 


1 Q A 


2 a 
Ar 


1.25 


3.16 


2.80 


0.98 


DoP 


10 


10 


10 


9 


Additional Results 


Initial R V i r (kpc) 


Final R V i r (kpc) 


Initial c v ir 


b mal c V ir 


NFW Halo 


150.0 


137.7 


10.2 


10.4 


Einasto Halo 


133.8 


175.0 


10.5 


9.9 


Reduced \ 2 an d degrees 


of freedom. 








'Final Einasto index n=3.1. 








Table 3: Rotation 


curve decomposition with 






M d =3.9 xlO 10 M Q (1 


ess 30% disk mass). 








Principal Results 


Hernquist 


Burkcrt 


NFW 


Einasto b 


Halo Mass (M©) 


8.5xlO n 


8.2xlO n 


1.6xl0 12 


2.1xl0 12 


Halo Radius (kpc) 


13.7 


3.0 


9.0 


16.8 


2a 
Ar 


0.98 


8.28 


0.96 


1.27 


DoP 


10 


10 


10 


9 


Additional Results 


Initial R V i r (kpc) 


Final R V i r (kpc) 


Initial c V i r 


Final ew 


NFW Halo 


150.0 


196.0 


10.2 


9.6 


Einasto Halo 


133.8 


215.6 


10.5 


9.4 



"Reduced \ 2 an d degrees of freedom. 
''Final Einasto index n=3.7. 



15 



Table 4: Maximum disk solution for NGC 5278 
M d =2.1 xlO 11 M Q . 
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surmount the RC maximum, practically the stellar 
disk does not fit any point of the observed RC of 
NGC 5278, so w e decided to use the definition of 



McGaughl ([20051 ) to explore the maximum disk so- 



lution for NGC 5278, increasing the baryonic disk 
we have determined in section[4]by 195% to fit the 
majority of the internal points of the observed RC 
of NGC 5278 without overtaking the observed RC 
maximum. We subtracted the determined maxi- 
mum disk from the observed RC of NGC 5278 to 
obtain the DM RC in the case of the maximum 
disk hypothesis and used the same DM halos fit- 
ting functions we have already used throughout 
this analysis to fit the DM RC. For this purpose 
we set the DM halos radius far from (~ 60 kpc) 
the last measured point (6.77 kpc) of the observed 
RC of NGC 5278. The results are shown in fig- 
ure [13] , [14] and the DM halos masses are listed in 
table[4] From the plots and the tabulated DM halo 
masses we can see clearly that there is still miss- 
ing a certain amount of mass (~ 30%) to repro- 
duce the entire mass of NGC 5278 enclosed in the 
observed RC. Classically, the ma ximum disk hy- 
pothesis ( van Albada et al.1ll985l ) have been con- 
sidered to explore the possibility to represent the 
RC of some spiral galaxies only with the contri- 
bution of the baryonic component. In the case of 
NGC 5278 we already know from section[4]that the 
baryonic component cannot account for the total 
mass of this galaxy even in the case of consider- 
ing more or less 30% errors in the determination 
of the baryonic mass. The DM analysis we have 
performed in section [5] tells us clearly that in the 
case of NGC 5278 a large amount of DM (» 3.75 
times the baryonic mass) component is required 
to entirely explain the observed RC of this galaxy 
and the maximum disk hypothesis confirms the 
conclusions of the RC decomposition study. 



6. Maximum Disk solution for NGC 5278. 



In the literature there does not exist a gen- 
eral definition of maximum disk. According to 



McGaughl (|2005l ) the maximum d isk M/L ratio 
should not exceed the RC data, and Gentile! ( 20081 ) 
defines the maximum disk as a fit where the stellar 
disk can totally account for the RC data peak. In 
the case of NGC 5278 the RC maximum is 342 km 
s _1 and if we increase the stellar disk in order to 
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Fig. 13.— Maximum disk solution with M d = 2.1 xlO 11 M Q . 
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7. Discussion. 

At the present epoch, interactions occur in a 
great variety of forms and environments. In this 
analysis we are interested in first stage interac- 
tions in order to deal with the evolution of baryons 
and DM during the interaction process of NGC 
5278 and its companion NGC 5279. In partic- 
ular, to study the evolution of baryons and DM 
halos during an interaction process is equivalent 
to first analyze the effects of the DM halo physics 
on the galaxy disks of the individual pair compo- 
nents and then to consider the consequences in- 
duced by tidal interactions. For this reason we 
begin dealing with the repercussions of the DM 
distribution on the formation and evolution of iso- 
lated galaxy disks and then we examine the whole 
scenario, adding the influences of tidal forces in 
shaping the DM distribution i n interacting galaxy 
pairs. The innovative work of Ostriker fc Peebles! 
(|l973l ) demonstrates that stiff spherical DM ha- 
los decrease the self-gravity of disks and stabi- 
lize them against gravitational instability (i.e. bar 
formatio n). On the o t her h and the more recent 
study of lAthanassoulal (|2002l ) establishes that the 
formation of strong bars, due to angular momen- 
tum transfer from the stellar disk, is favored when 
the DM halo is simulate d with collisionle s s par - 
ticles and axi-symmetric. 



Berentzen etaD (|200d ) 



proved that triaxial DM halos create chaos in a bar 
, reducing the bar instability. DM halos can effi- 
ciently slow dow n bars through dynamical friction 
(jWeinbergl Il985j ) and some authors infer the DM 
conten t of barred galaxies from the bar p attern 
speed (|Debattista fc Sellwoodl Il998l l2000h . The 
misalignment of the baryonic infall with the minor 
axis of the DM halo removes a large amount of an- 
gular momentum from baryons during galaxy for- 
mation; as a consequence only co mpact spheroids 
can form (|Aumer fc Whitell2012l ). 
From the first self-consistent numerical simula- 
tions of disk/halo systems it was became clear 
that the DM halo play a fundamental role to ex- 
haust the binding energy and angular momen- 
tum of the luminous components changi ng dra- 
matic ally its distribution and dynamics ( Barnes! 
1988). Successive numerical studies of galaxy 
interactions attempt to constrain the DM halo 
mass distribution to the morphology of the tidal 
tails that generate during t he dynamical evolutio n 
of the interaction process ( Dubinski et al. 1996). 



Mihos et al.l (|1998l ) explore the aftermath of dif- 
ferent DM halo potentials on the kinematics and 
morphology of tidal tails considering a dynamical 
model of NGC 7252 and find that halo masses of 
the order of 4-8 times the baryonic mass compo- 
nent (disk+bulge) reproduce the observed interac- 
tion scenario of NGC 7252; nevertheless, the same 
authors note that numerical simulations with cos- 
mological initial conditions are needed to entirely 
validate the use of tidal tails to constrain the mass 
of D M halos in the evolut i on of the interaction pro- 
cess. ISpringel fc White! (|l999j ) perform the tidal 



tails test in the framework of CDM models and 
obtain that in the context of the analyzed CDM 
models, it is difficult, that tidal tails could con- 
strain DM halo parameters. ISpringel fc Hernquist 
(|2005f ) study the remnants of interaction processes 
in the final stage of evolution at high redshift 
of gas-dominated galaxies by means of numeri- 
cal simulations considering a disk and a DM halo 
without any bulge to determine that the rapidly 
cooling of the gas feeds a star-forming disk resem- 
bling a spiral galaxy and not an early- type object. 
The subject of the evolution of DM halos in inter- 
actions is too vast to allow a detailed treatment of 
the theme in this discussion, we have only briefly 
addressed here the most relevant points that are 
of more interest for the analysis we have made 
in this work. In particular we note that several 
results are still under debate and many aspects of 
the evolutionary problem of the DM halo in galaxy 
interactions are at the time totally unknown. 
In this study we offer a very modest perspective 
that analyzes the interplay between the baryonic 
component of NGC 5278 and the DM halo of 
this interacting galaxy (only at the present time) 
through the RC decomposition technique, consist- 
ing in fitting to the galaxy gravitational potential 
different DM components, once the baryonic pa- 
rameters are determin e d wit hin a robust method 
(e.g., the lZibetti et al.l ([20091 ) method). 

We perform here a simple and approximate es- 
timate of the tidal radius and dynamical time for 
KPG 390 in order to assess the validity of our 
RC analysis. In our former work we have already 
determined the interaction time for KPG 390 to 
be T int = 2.9 x 10 s yrs. We consider NGC 5279 
as a spherical mass distribu tion with total mass 
m s = 2.2 x 10 10 according to lLequeuxl |l983), and 
NGC 5278 is composed by an exponential disk and 
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Fig. 15. — Plot of the dynamical time vs tidal radius for KPG 390. The dashed line indicates the interaction 
time for KPG 390. 



a DM halo with mass M{R) = M d (R) + M h (R), 
where M d (R) is the exponential disk mass, Mh{R) 
represents the DM halo mass, and R kpc is our RC 
extension. We consider the best fitting halos of our 
analysis, the Hernquist DM halo and the Einasto 
DM halo. The tidal radius is defined accord- 
ing to the relation R t = Xi 2 (m s /(fcM(i?))) 1/3 , 
where X\2 — 16.6 kpc is the separation betwe en 
the two pair component (|Karachentsevl fjl972h ). 
and k — 4 indicates a circular orbit. The dy- 
namical time is defined according to the relation 
t d = (2R z t IGM(R t )) 1 / 2 , where R t represents the 
tidal radius. The results are shown in figure 1151 
and we can see that the minimum tidal radius 
begins around 6.8 kpc in the case of the cored 
(n < 4) Einasto DM halo and more 30% disk 
mass, while the observed RC of NGC 5278 ex- 
tends up to 6.7 kpc; thus in principle the tidal 
perturbation induced by the companion of NGC 
5278 does not affect so much the kinematics of 
NGC 5278 described by the Ha RC. The dynam- 
ical time crosses the interaction time line beyond 
20 kpc, in the case of Hernquist DM halo with 
disk mass M d = 5.6 x 10 10 M Q and less 30% disk 
mass, so that the interacting system enters the 
zone of interaction far beyond the extension of the 
Ha RC. From this rough guess we conclude that 
the RC decomposition performed in this study is 



not so deeply influenced by the ongoing interac- 
tion process between NGC 5278 and NGC 5279. 
The observed Ha RC of NGC 5278 reaches 6.8 kpc 
(« 13 arcsecs) and the baryonic disk extends up to 
17.3 kpc (w 33 arcsecs). The RC decomposition 
analysis gives us some prescriptions about which 
kind of DM halos better fit the observed RC of 
NGC 5278, nevertheless the short span in radius 
(rs 1/2.5 of the baryonic disk radius) of the RC 
of NGC 5278 forces us to consider these results 
only as a first approximation. In fact we ignore 
which could be the behavior of the RC of NGC 
5278 beyond the last measured point of our Ha 
RC. On the other hand, it is interesting to note 
that in the case of disk mass M d — 5.6 x 10 10 M@, 
if we consider the Hernquist best fitting DM halo 
with final mass M d = 8.7 x 10 11 M Q , we need 3.75 
times the mass of the baryonic disk to account for 
the entire mass of NGC 5278. The necessity of 
such a large amount of DM for NGC 5278 within 
a radius of only 6.8 kpc is surprisi ng given that for 
instance Ivan Albada et al.l (|1985l ) found a similar 
DM amount (4 times the disk mass) in NGC 3198, 
but within 30 kpc. The same authors observe a 
decrease of the amount of DM (1.5 times the disk 
mass) inside 15.9 kpc. If we extrapolate that ob- 
served non linear trend to the RC of NGC 5278, we 
could expect much more DM if the observed RC 
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could reach the baryonic disk extension, assuming 
a flat RC up to the optical outskirts of NGC 5278. 
From these considerations we can state that even 
if the short range of the Ha RC of NGC 5278 does 
not allow us to derive general properties of DM, 
notwithstanding it reveals the large DM content 
of NGC 5278 within just a few kpc. 
In the Introduction we have reminded that 
KPG 390 is classified bv lKlimanov fe Reshetnikov 
(|200l[ ) as M51 galaxy type. In the literature 
there are few examples of numerical modeling 
of these systems including the archetype M51. 
In particular, we are interested in trying to fit 
the coupled evolution of DM and baryons dur- 
ing at least the present epoch of the interaction 
process of M51-type interactions. For this rea- 
son we concisely recapitulate here the most im- 
portant steps, beginning with the first attempts 
to reproduce this interacting systems, up to the 
actual state-of-the-art numerical experiments of 
this kind of interactions. The landmark study 



of iToomre fc Toomrd ([19721 ) succeeded to repro- 
duce the principal morphological features of M51, 
including bridge and tails, considering parabolic 
encounters of two galaxies idealized as disks of 
noninteracting tests particles without self-gravity. 
This pioneering work was unable to reproduce the 
grand-design spiral structure of M51 and this left 
out was a t tribute d to the despise of self-gravity. 
Hernquist (1990b) was the first to accomplish self- 
consistent 3D numerical simulations of M51 and 
focused attention on parabolic orbits. The au- 
thor pointed out that the incorporation of self- 
gravity is not enoug h to gener a te th e extended 
Hi-tail discovered bv lRots et al. ( 199Cll ). instead a 
longer duration of the perturbation is neces s ary to 



recreate the elongated HI feature. iToomrd ((1994) 
improved the original model of M 51 adding self- 
gravity and increasing the duration of the satel- 
lite crossing since the perturbation. The result 
was a better morphological fea tures m a tching , no- 
tably for the large HI- tail. iBarnesI (1998) has 
noticed that the interpretation of the M51 inter- 
action is far from being adequate, principally for 
the seemingly counter-rotating direction of the 
extended Hl-tail, that probably denotes signifi- 
cant inclination with respect to the inn er disk 
of M51. For this reason, according to Barnes! 



Salo fc Laurikainen (1993) numerically simulate 
the M51-type pair ARP 86 (NGC 7753-7752) con- 
sidering each galaxy composed of self-gravitating 
stars+gas disks embedded in rigid spherical ana- 
lytical DM halos. The resulting models match the 
morphology and kinematics of the principal tidal 
features and return right projected separations 
and radial velocity difference between the two pair 
co mponents. A major draw back of the approach 
of lSalo fc Laurikainen! (|l993l ) is that the DM halos 
are treated as rigid components, thus the transfer- 
ence of material between each component of the 
pair cannot produce any effect on the DM distri- 
bution of the system, and on the other hand the 
DM distribution does not vary during the inter- 
action process, fact that may affect the baryonic 
infall of one component above the other. These 
authors apply the same tech niques in a more re- 
cent n umerical study of M51 ( Salo fc Laurikainen! 
2000al lbT). They succeeded to reproduce the outer 
(Hl-tail) and inner morphological features of M51 
and in the case of these studies the choice of rigid 
DM halos is not so crucial because the authors 
are not specifically interested in a detailed study 
of the relation between baryons and DM dur- 
ing the interaction process, yet this methodology 
leaves open many questions about the actual in- 
terplay b etween the m a in ma ss components in 
galaxies. Rosado et al. ( 201lh perform numeri- 
cal simulations of the M51 galaxy type KPG 302 
(NGC 3893/96), where NGC 3893 is a grand- 
design Sc type galaxy and N GC 3896 is cataloged 
as a l enticular galaxy (SOa). IVerheiien fc Sancisl 
(|200ll ) by means of HI imaging demonstrated 
the existence of a prolonged envelope surround- 
ing the two galaxies of the pair in the direction 
SE to NW alongside the imaginary line connect- 
ing the nuclei of both components of KPG 302. 
KPG 390 is morphologically somewhat different 
from KPG 302, nevertheless we briefly summa- 
riz e here the results of the numerical simulation 
of Rosado et al. I (l20ll to acquire more insight 
about the dynamical evolution of KPG 390. The 
numerical simulati ons of KPG 302 were achieved 
with GADGET2 dSpringel et all l200lt ISnringell 



2005) and using a galaxy mod el analogous to that 
of iBarnes fc Hernquist! ( 19961 ). In these models 



(|1998l ) , it is complicated to produce this tilt in the 
framework of the above-quoted numerical models. 



stars and DM are treated as collisionless fluids, 
the gas is regarded as a compressible isothermal 
fluid that undergoes gravitational forces and hy- 
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drodynamical forces as well. The gas forms 10% 
of the stellar disk and shocks are managed in- 
cluding an artificial viscosity in the equations of 
motion. The gaseous and stellar disks have expo- 
nential density profiles and the DM density profiles 
are represented by Hernquist, P-ISO and NFW. 
The satellite galaxy NGC 3896, is represented by 
means of a Plummer sphere. The orbits consid- 
ered for the simulations are parabolic. The results 
of the simulations resemble the observed morphol- 
ogy of KPG 302 and the measure d Ha velocity 
maps (|Fuentes-Carrera et al.l 120071 ). The authors 
outlined, it seems from the simulations that with 
P-ISO DM halos more open spiral patterns are 
readily formed. This fact, if confirmed with more 
numerical tests, probably indicating strong dif- 
ferences in shaping spiral patterns, depending on 
the DM density profile considered. This occur- 
rence could be connected with the deepness of the 
DM halo potential well, given that cored DM den- 
sity profiles have shallower potential wells than 
cuspy DM profiles; the matter in the spiral arms 
is rearranged at different azimuthal distances in 
the galaxy disk depending on the gravitational 
response of the spiral arms baryons to the col- 
lisionless DM particles. If this hypothesis holds 
at least for M51-type interactions, we would ex- 
pect for NGC 5278 a less open spiral pattern due 
to the cuspy Hernquist (in the case of disk mass 
M d = 5.6 x 10 10 M Q and less 30% disk mass) DM 
halo gravitational potential well, i. e. we could 
constrain a range for the pitch angle of the spiral 
arms of this galaxy through the knowledge of the 
DM distribution. In the case of more 30% disk 
mass the success in the RC fitting of the cored 
(n < 4) Einasto DM halo depicts a very different 
story in which probably the spiral pattern of NGC 
5278 has a wider pitch angle with respect to the 
Hernquist DM halo solution. At the present time 
these latter considerations are only on a conjec- 
tural ground and it is clear that a much more work 
is needed to deeply understand the mutual reci- 
procity between DM and baryons in interacting 
spiral galaxies. 

In particular, in the case of KPG 390, it is difficult, 
at the present moment, without a numerical sim- 
ulation of this pair, to establish an evolutionary 
path for the DM halo distribution in this system, 
due to the uncertainties about the orbital configu- 
ration parameters that determine the geometrical 



evolution of this interaction. Successive passage 
of the companion NGC 5279 could diminish the 
baryonic content of NGC 5278 leading to very dif- 
ferent configurations and shapes of the DM com- 
ponent in this galaxy. Summarizing, we can affirm 
that the present analysis does not exhaust the dy- 
namical study of KPG 390, but settles important 
prescriptions about the DM content of NGC 5278, 
the primary galaxy of KPG 390. We plan to use 
this information as boundary conditions in future 
numerical simulations of this complex interaction 
process. 

8. Conclusions. 

In this article we revisit the mass distribu- 
tion of the spiral galaxy NGC 5278, this time 
with new observational information about the disk 
M/L coming from surface photometry and stel- 
lar population synthesis studies. We derive 2D 
M/L stellar maps and find that the M/L varies 
along the galactocentric radius. We perform ID 
and 2D bulge-disk decomposition to conclude that 
our analysis does not reveal any bulge component 
in NGC 5278, alternatively we find out a powerful 
central source and a classical exponential disk. We 
obtain the 2D disk mass distribution of NGC 5278 
and from this distribution we derive the disk mass 
contribution to the whole RC of NGC 5278, con- 
sidering a disk mass variation of mor e or less 30%, 
to inc lude the uncertainties of the IZibetti et al 



( 2009t ) method. We realize that the baryonic RC 
of NGC 5278, and also the baryonic RC more and 
less 30% disk mass, is not able to explain the whole 
mass content of NGC 5278 revealed from its gravi- 
tational potential. For this reason we subtract the 
baryonic disk of NGC 5278 from the observed RC 
of this galaxy, also in the cases of more and less 
30% disk mass, to obtain the DM RCs for NGC 
5278. As a next step we fit to the DM RCs of 
NGC 5278 four different DM halo distributions to 
establish which of the four DM halos better repro- 
duce the DM RCs of KPG 390 A. The main results 
of this work are that the Hernquist DM halo bet- 
ter fit the DM RC of NGC 5278 in the case of 
disk mass M d = 5.6 x 10 10 M Q and less 30% of 
this value, and also that the cored (n < 4) Einasto 
DM halo better fit the DM RC of NGC 5278 in the 
case of more 30% disk mass. We also try the maxi- 
mum disk solution for NGC 5278, fitting the maxi- 
mum number of internal points of the observed RC 
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with 195% of the determined baryonic disk. The 
maximum disk solution is inadequate to account 
for the entire mass of NGC 5278. We are unable 
to predict the time evolution of the baryonic and 
DM distribution because we are missing a whole 
dynamical study of the pair KPG 390. For this 
reason we plan to use these results on the mass 
distribution of NGC 5278 as boundary conditions 
in a future numerical simulation of the pair KPG 
390 to unveil other important dynamic features of 
this complex type of interaction. 
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